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We propose a method for the theoretical investigation of polymer transloca- 
tion through composite pore structures possessing arbitrarily specified geometries. 



Translocation through each constituent part of the composite is treated as being 
analogous to the diffusion of the translocation coordinate over the free energy land- 
scape derived from the chain configurations within the pore. The proposed method 
accounts for possible reverse motions of the leading chain end at the interface between 



constituent parts of a composite pore, a possibility that has been neglected in prior 
studies. As an illustration of our method, we study the translocation of a Gaussian 
chain between two spherical compartments connected by a cylindrical pore, and by 
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00 ' a composite pore consisting of two connected cylinders of different diameters, which 

o ' 

is structurally similar to the a-hemolysin membrane channel. We demonstrate that 
reverse chain motions between the pore constituents may contribute significantly to 
the total translocation time. Our results further establish that translocation through 
a two-cylinder composite pore is faster when the chain is introduced into the pore 
on the cis (wide) side of the channel rather than the trans (narrow) side. 
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I. INTRODUCTION 

The migration of biopolymers through nanopores plays a key role in several biological 
processes, including the transport of mRNA molecules from the nucleus to the cytoplasm 
following transcription, the transport of proteins to and from the nucleus, and the injection 
of viral DNA into a host cell 1 . Biopolymer translocation further has technological appli- 
cations in the development of biosensors for polynucleotide analysis and sequencing^^^. 
In view of its biological and technological significance, polymer translocation has been the 
subject of several recent experiment a^yjlkLSAifi theoretical^^^^^^^^^ 1 ^ 1 ^ and 
simulation-based^ 1 ^ 1 ^ 1 ^ 1 ^ 7 - 1 ^ 1 ^ investigations. These, as well as related studies concerning 
polymer translocation phenomena, are summarized in the review of Meller— . 

While many prior theoretical investigations have focused on translocation through a nar- 
row pore embedded in a rigid wall of negligible thickness, wherein the pore is assumed to 
accommodate at most a single polymer segmentii^, pores in biological systems can be 
significantly large. Recently, there has been much interest in examining the effect of pore 
geometry on the translocation process. Studies in this direction include investigations into 
the influence of pore length on the translocation time^ and the release of a chain from 
one spherical vesicle into another—. Pore structure is expected to play a significant role in 
controlling translocation through membrane channels in biological systems. In particular, 
experimental studies have shown that translocation occurs more readily when the chain is 
introduced through the side of the asymmetric a-hemolysin membrane channel that has the 
larger vestibule 7 . An investigation of the influence of pore geometry on translocation dynam- 
ics entails the detailed consideration of chain configurations within the pore, as well as chain 
motions through the constituents of a composite pore structure. To this end, Muthukumar— 
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has considered polymer translocation through several pore geometries, including the migra- 
tion of a chain from one sphere to another through a narrow cylindrical pore. The latter 
study was further generalized by Wong and Muthukumar 22 , who took account of the finite 
diameter of the cylindrical pore by allowing the number of segments within the pore to vary 
and explicitly determining the corresponding free energies of the confined chain. 

The works of Muthukumar— and Wong and Muthukumar— provide a framework for the 
explicit treatment of polymer translocation through composite pore structures with account 
for all possible chain configurations within the pore. However, the assumptions underlying 
these studies do not allow for the unrestricted motion of the chain at the interface between 
adjacent constituent parts of the pore. Specifically, during translocation from the donor 
sphere to the receptor sphere through a cylindrical pore, the leading chain end was assumed 
to be "absorbed" at the cylinder-receptor sphere interface immediately upon its arrival at the 
cylinder end, and reverse chain motions at the interface were disallowed 17 . A subsequent 
attempt was made to relax this assumption, by allowing a small, fixed number of chain 
segments at the leading end to retreat into the cylinder even after entering the receptor 
sphere 22 . Nonetheless, possible chain motions at the interface continue to be curtailed under 
this assumption, whereas, in reality, the occurrence of reverse chain motions may contribute 
significantly to the total translocation time. 

In the present contribution, we propose a modification of the method of Muthukumar and 
coworkers^ 7 - 1 ^ 2 - whereby the leading chain end is allowed to move backward as well as forward 
at the interface between adjacent parts of the pore, and may arrive at the interface any num- 
ber of times before chain migration into the pore constituent immediately downstream is 
initiated. This is achieved by imposing the more realistic radiation boundary condition^ in 
place of the absorbing boundary condition on the translocation coordinate at the interfaces 
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between constituent parts of the pore. In the case of chain transport from one sphere to 
another through a cylindrical pore, we demonstrate that backward chain motions may lead 
to a significant increase in total translocation time. Our results indicate that the transloca- 
tion time increases monotonically with increase in cylinder length, rather than exhibiting a 
minimum with respect to the cylinder length as predicted by Wong and Muthukumar—. We 
next consider the translocation of a chain from one sphere to another through a composite 
two-cylinder pore, which serves as a simple model of the mushroom-shaped a-hemolysin 
membrane channel. Our results confirm that translocation is faster when the chain is in- 
troduced into the pore on the cis (wide) side of the channel rather than the trans (narrow) 
side. 

The paper is organized as follows. In Sec. HIl we treat the translocation of a Gaussian 
chain from the donor sphere to the receptor sphere through a cylindrical pore. Translocation 
through an a-hemolysin pore is considered in Sec. Ifflt while Sec. IIVI provides a summary 
of our findings. 



We first briefly summarize the approach previously employed by Park and Sung and 
Muthukumar and coworkers^ 1 ^ and adopted here by us for the derivation of the free energy 
landscape for the translocation process. The probability density P(r, r , N) that the ends 
of a Gaussian chain of N segments are at positions r and r is governed by the equation^ 



II. 



TWO SPHERES CONNECTED BY A CYLINDER 



A. Free energy landscape 




(1) 



where the unit of length has been set equal to the Kuhn length (denoted b by Wong and 
Muthukumar). Equation (JTj) is solved subject to the condition P(r, r , 0) = 5 3 (r — r ) and 
the condition that P(r, ro, N) vanish at all surfaces. 

The solution of Eq. (TjQ) in conjunction with the appropriate boundary conditions for a 
chain confined to a sphere of radius R yields the expression 
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with the origin chosen to lie at the sphere center. If one chain end, say, r , is allowed 



to lie anywhere within the sphere while the other end, r, is tethered arbitrarily close to 
the sphere surface such that r = R — c, with c the tethering distance, then the probability 
density per unit area of the sphere surface corresponding to r becomes 
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Similarly, we obtain the expression 
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for a chain confined to a cylinder of length M and radius a, and whose ends lie at positions 
r = (r, 9, z) and r = (r , 6* , z ) measured relative to an origin located at the center of the 
cylinder axis. If one chain end, r , is free and the end at r is fixed such that z = —M/2 + c, 
where the z-coordinate is measured along the cylinder axis and c is a small tethering distance, 
then 
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For a chain whose ends are fixed such that z = —M/2 + c\ and z = M/2 — c 2 , we obtain 
the expression 

87r 3 a 2 cic 2 ^ m 2 +1 [ f f4 >k , m 2 n 2 \ N 

The above expressions are correct to leading order in the tethering distances, and the symbols 
^o,k, k = 1,2, ... appearing therein are the roots of the Bessel function of the first kind of 
order 0, denoted by Jo- We note that the numerical prefactors in Eqs. (JI])-© differ slightly 
from those of Wong and Muthukumar—. 

Equations ([3]), ([5]) and ([6]), when appropriately combined, enable the computation of the 
free energy landscape for the translocation of a chain from the donor to the receptor sphere 
through a cylindrical pore, as illustrated by the sequence of chain configurations in Fig. CD 
As evident from Fig. [TJ we restrict our attention to situations where the chain length greatly 
exceeds the pore length, and disregard situations wherein the chain completely enters the 
cylinder and subsequently diffuses to the receptor sphere entrance. We will elaborate on this 
assumption in the following section. Following Wong and Muthukumar, we denote by p and 
q the number of chain segments in the cylinder and in the receptor sphere, respectively. We 
further denote the number of segments in the cylinder just before the start of translocation 
into the receptor sphere [cf. Fig[T](c) and (d)] by R p , with M < R p < N. The radii of the 
donor and the receptor spheres are denoted, respectively, by R± and i?2- The free energies 
of the chain configurations depicted in Fig. [T]are22 (in units of fcsT) 
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F (a) = -ln [P Sl (Ri,N)7ra 2 ] 
F {b) (p <R P ) = - In [P Sl (Ri, N - p)P Cl (p)] 
F {c) (p = R p <N) = -\n [P Sl (R u N - p)P c . 2 (p)} 



(7) 
(8) 
(9) 



F (d) (p = R p = N) = - In [P C2 (N)\ (10) 

F (e) (g, R p <N-l) = -\n [P Sl (R u N — Rp — q)P C2 (R p ) P Sl (R 2 , g)] (11) 

F (f) (q, R, < N) = - In [P Ca - g) P 5l (i? 2 , g)] (12) 

F W (g) = - In [P Cl (N - q) P Sl (R 2 , g)] (13) 

F {h) = -ln[P Sl (R 2 ,N)7ca 2 } (14) 

In the subsequent analysis, all tethering distances are set to 1/2 (in units of b) for 
simplicity 17 . 

B. Translocation time 

The sequence of steps constituting the translocation process may be divided into two 
stages for the sake of convenience. The first stage, which begins with the insertion of the 
first segment of the chain into the cylinder from the donor sphere and lasts until the initiation 
of translocation into the receptor sphere, is comprised of either the sequence (a)— »-(b)— >(c) 
for R p < N or the sequence (a)— ►(b)— »(d) for R p = N in Fig. [TJ The transport of the chain 
into the receptor chamber occurs entirely within the second stage of translocation, which 
may be comprised of either the sequence (c)— ►(e)— >(f)— >(g)— >(h) for R p < N [where step 
(e) must be omitted if R p = N — 1] or the alternate sequence (d)— >(g)— >(h) for R p = N. 

The probability density W p (t) of the number of segments p contained within the cylinder 
at time t during the first stage of translocation is governed by the equation 34 



dW p _ d_ 
dt dp 



(15) 



dp p dp 

where we have selected the unit of time to be the reciprocal of the diffusion constant ko (in 
the notation of Wong and Muthukumar), and F denotes the free energy of the chain [cf. 
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Eqs. (I7|)- (|T4|) ]. Consistent with prior theoretical studies, we decouple the diffusion of the 
chain end from within the donor sphere to the cylinder entrance from the process of chain 
migration into the cylinder once the leading chain end has located the pore. Therefore, we 
impose a reflecting (no flux) boundary condition at p = 1: 



OF 9W P 



= (16) 

j p =i 



dp p dp 

The number of segments R p in the cylinder at the end of the first stage of translocation 
must be at least M and is at most N [cf. Fig[T](c) and (d)]. Moreover, the leading end of the 
chain may arrive at the cylinder-receptor sphere interface many times before finally being 
"absorbed," thereby initiating the second stage of translocation. For this reason, we impose 
the following radiation boundary condition 3 -^ in place of the absorbing boundary condition 
at p = R p : 

k R W Rp (t) (17) 



dp Wp+ dp Jp=Rp 

where k R is a "reaction" rate constant, expressed in units of k . The radiation boundary 
condition reduces to the absorbing boundary condition in the limit k R —>■ oo, and to the 
reflecting boundary condition in the opposite limit k R — > 0. The average passage time 
to reach p = R p starting from p = 1 for a process governed by Eq. (Tl5|) and boundary 
conditions (fT6l) and (TTT1) is given by the expressior 31 



n(R P ) = J " dy'exp [F(p = y')] J" dy" exp [-F(p = y")} 



+ -L exp [F(p = Rp)} dy'exp [-F(p = y')\ (18) 



We emphasize that the chain configurations illustrated in Fig. [T]are exhaustive only under 
the assumption iV 3> M. We enforce this condition so as to avoid situations wherein the 
chain completely enters the pore and subsequently diffuses across the length of the pore to 
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the receptor sphere entrance. The latter situation, which arises when N <C M, necessitates 
the augmentation of the translocation time by the inclusion of a diffusion time in addition 
to the time taken by the chain to completely cross into the cylinder. We here restrict our 
analysis to the physically realistic situation wherein N S> M. 

Likewise, the probability density W q (t) of the number of segments q present in the receptor 
sphere at time t during the second stage of translocation is governed by the equation 



dW q _ d_ 
dt dq 



sf aw. 



(19) 



dq ' q dq 

subject to a reflecting boundary condition at q = 1 and an absorbing boundary condition 
at q = N. The latter absorbing boundary condition is motivated by the fact that only 
successful translocation events are recorded, while the reflecting boundary condition at q = 1 
is consistent with the use of a radiation boundary condition at p — R p , signifying the 
initiation of translocation into the receptor sphere. The corresponding mean first passage 
time is^ 

/N—l ry' 
dy'exv[F(q = y')]J dy" exp [-F(q = y")\ (20) 

For given R p , Eqs. (ITS1) and (120]) yield the corresponding times of translocation, Ti(R p ) 
and T2(Rp). We equate the average durations (ti) and (r 2 ) of the first and second stages 
of translocation, respectively, to the average of Tx(R p ) and T2(R P ), respectively, computed 
with respect to the probability of configurations having R p segments in the cylinder for 
Rp = M, M + 1, ...N [cf. FigG^c) and (d)] as follows: 



,. , _ Er;1mP Si (Ri,n - R p )P C2 (R P )n(R P ) + Pc 2 (N)n(R p = N) i 
J2t=M P s 1 (Ri,N-R P )Pc 2 (R P )+Pc 2 (N) 



and 

yN-l 



, _ x _ EZ=m P sARi,N- R p )P C2 (R p )t 2 (R p ) + P C2 (N)t 2 (R p = N) ^ 
""' 2) Y!1:1mPsARi,N-R p )Pc 2 {R p ) + Pc 2 {N) 



10 



The total time of translocation may now be defined as (r) = (ti) + (T2). 

The rate constants which are dependent on the value of R p , remain to be specified. 
These are derived from the radiation boundary condition given by Eq. f TTTT) in conjunction 
with the backward finite difference approximation 

- In [P Sl (R u N - R P )P C2 {Rp)] + In [P Sl {R 1 ,N-R p + l)P Cl (R P - 1)] , R p < N 

- In [P C2 (N)} + In [P Sl (Ru l)P(h(N - 1)] , R p = N 

(23) 



dF 
dp 



P=R P 



We further replace the term W p (t) appearing in Eq. (TTTT) with the probabilities of the 
corresponding configurations: 

P Sl (Ri,N -R P )P C2 (R P ), R P <N 



W p (p = R p ) ~ < 

Pc 2 (N), R P = N 

V 

and the backward finite difference approximation yields 



(24) 



dW v 



dp 



P=R P 



P Sl {R u N - R P )P C2 {R P ) - P Sl {Ri, N-R p + l)P Cl (R P ~ 1), R P < N 
PaW-P^RuVPaiN-l), R P = N 



(25) 



We note that the simultaneous derivation of F and W p from the probabilities of the chain 
configurations would result in a vanishing k R for a continuous process [cf. Eq. fTTTT) ] . How- 
ever, in the present case of discrete chain transport in the absence of external forces or fields, 
the diffusive contribution to the probability flux in Eq. (TTTT) dominates over the convective 
contribution, owing to the large change in the probabilities of the corresponding chain con- 
figurations when p changes from R p — 1 to R p . In our calculations, values of Jsr obtained 
were in the range of O(l)-O(100). 

Use of Eqs. (121]) and fl22]) . following the numerical integration of ffl8]) and fl20l) in combi- 
nation with (jTf)- (jT4"]) . fTTTT) and fl23]) - fl25]) . finally yields the desired translocation times. 
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C. Results 

Our results for a chain oi N = 300 segments are illustrated in Figs. EHU where length 
and time have been expressed in units of b and A;^ 1 , respectively. Figures [2] and [3] depict 
translocation times as a function of M for several values of R 2 with Ri = 30, and a = 3 
and a = 5, respectively. It is evident that, contrary to the observations of Wong and 
Muthukumar—, the total translocation time grows monotonically with respect to increasing 
cylinder length. The results suggest that backward motions at the cylinder-receptor sphere 
interface do, in fact, significantly contribute to the duration of the first stage of translocation. 
The chain incurs an entropic penalty upon moving into the cylinder from the donor sphere, 
whereas there is an increase in entropy during translocation into the receptor sphere. Hence, 
the total translocation time (r) is dominated by the contribution from the first stage of 
translocation and, consequently, reflects the increase in (n) with increasing M. 

The duration of the second stage of translocation decreases as M is increased, owing to 
the fact that chain configurations in which a larger number of segments are present in the 
cylinder have a higher free energy than those in which a larger number of segments are 
present in the donor sphere. As a result, translocation into the receptor sphere leads to a 
greater free energy drop as M is increased. This observation was made earlier by Wong and 
Muthukumar. However, the second stage of translocation makes a far smaller contribution 
to the total translocation time relative to the first stage. 

Figure H] illustrates the translocation times as a function of M for several values of a, 
with R\ = 30 and R2 = 60. As expected, an increase in a causes a decrease in (r) and (ti), 
owing to the concomitant lowering of the entropic barrier to translocation into the cylinder. 
Moreover, configurations in which a larger number of chain segments reside in the cylinder 
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become more favorable as a is increased relative to those in which the majority of chain 
segments lie in the donor sphere. Since the former configurations have a higher free energy 
than the latter, the driving force for translocation into the receptor sphere is greater in the 
former case. Consequently, fa) decreases with increase in a. 

III. TWO SPHERES CONNECTED BY A COMPOSITE TWO-CYLINDER 

PORE 

In this section, we consider the translocation of a chain from a donor sphere of radius R\ 
to a receptor sphere of radius R2, connected by a composite two-cylinder pore with radii a% 
and ct2 and lengths M\ and M2, where the subscripts 1 and 2 refer to the outer and the inner 
cylinder, respectively. We now divide the translocation process into three stages, namely, (1) 
the transport of the chain into the outer cylinder until the leading chain end enters the inner 
cylinder, (2) the transport of the chain into the inner cylinder until its leading end enters 
the receptor sphere, and (3) the complete migration of the chain into the receptor sphere. 
We denote by I, m and n, the number of segments in the outer cylinder, inner cylinder and 
receptor sphere, respectively, and by Ri and R m the maximum number of chain segments in 
the outer cylinder and inner cylinder, respectively, just prior to the start of the succeeding 
stage. Again, we have the conditions that Mi < Ri < N, and M2 < R m < N. 

Under the assumption that the chain length greatly exceeds the lengths of the cylinders, 
Fig. [5] illustrates all possible chain configurations during translocation. The first stage 
of translocation comprises either the sequence (a)— >(b)— >(c) for Ri < N, or the sequence 
(a)^(b)^(d) for Ri = N (cf. Fig [5]). Similarly, the second stage of translocation is iden- 
tified with one of the following sequences: (c)— >(e)— >(f) for Ri + R m < N; (c)— >(e)— >(g) 
for R t + R m = N; ( c )->(e)->(h)->(i)->G) for R t < N, R m < N and Ri + R m > N; 
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(c) ->(e)->(h)->(i)->(k) for R t < N, R m = N; (d)->(i)->(j) for R t = N and R m < N; and 

(d) — >(k) for Ri — N and R m = N. Finally, the third stage of translocation is described 
by one of the following: (f)— >(m)— >(n)— >(o)— >(p)— >(q) for Ri + R m < N (where the 
intermediate step (1) must be omitted if Ri + R m = N — 1); (g)— >(n)— >(o)— >(p)— >(q) for 
Ri + R m = N; (j)->( n )->(o)->(p)->(q) for R m < N and R t + R m > N; and (k)->(p)->(q) 
for R m = N. The free energies of these chain configurations are listed in the Appendix. 

The probability densities Wi, W m and W n of the translocation coordinates /, m and n, 
respectively, are governed by the equation 



dW x 
dt 



d_ 

dx 



OR dW x 
— W x H 

dx dx 



(26) 



where x — I, m or n, and F denotes the free energy of the chain. We again impose a 
reflecting boundary condition at / = 1 and an absorbing boundary condition at n — N, 
corresponding to experimental measurements wherein only successful translocation events 
are recorded. Since the leading chain end may visit the interface between the inner cylinder 
and the outer cylinder or the outer cylinder and the receptor sphere many times before 
being "absorbed" at the boundary, we impose the radiation boundary condition at I — Ri 
and m = R m , given by the expressions 
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di 1 di 



hWi(l = R t ) 



(27) 



l=Ri 



and 



OF dW m 
om om 



k 2 W m (m = R m ) 



(28) 



m=R m 



where k\ and k 2 are rate constants (in units of k ), determined following the procedure 
described in the Appendix. In addition, we impose reflecting boundary conditions at m — 1 
and n — 1, consistent with the above radiation boundary conditions. 
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The corresponding average passage times are yielded by the expressions 

n(^) = j* dy'exp [F(l = 2/')] f dy" exp [-F(l = y")} 

+ 1 exp [F(Z = i*)] ^ ^'exp = </)] (29) 

/Rm ry' 
dy , exp[F(m = y')} J dy" exp [-F(m = y")} 

+ - exp [F{m = R m )\ / dy' exp [-F(m = y')\ (30) 

and 

/iV— 1 /-J/' 
dy'exp[F(n = y')]J^ dy" exp [-F(n = /)] (31) 

We finally obtain the average translocation times from the expressions 

^ = T,r 1=Mi Efl m =A/ 2 p (gt) P(R m \Ri)r 2 (Rl, R m ) ^ 

Y.R l =M 1 Y.R m =M 2 P ( R l) P ( R rn\Rl) 

and 

(r 3 ) = ^= M ^^=ivi2- ,-„.,-.,-., (34) 



Hr 1= Mi J2 Rm =M2 P W P ( R m\Rl) 



where P(Ri) and P(R m \Ri) denote the probabilities of configurations with a maximum of 
Ri segments in the outer cylinder, and a maximum of R m segments in the inner cylinder 
given the value of R[, respectively 

Equations (|27|) - (|34|) enable the computation of the translocation times for each of the 
three stages, whereby we obtain the total translocation time (r) from the sum (n) + (t 2 ) + 
(r 3 ). The details of the calculation are analogous to those presented in Sec. [TT1 and are 
relegated to the Appendix. 

Our results are illustrated in Fig. [6HH1 for a chain possessing N = 100 segments. Figure 
[61 which depicts the total translocation time for several values of a 2 with a% — 3 and a% = 4, 
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and vice versa, and with Mi = M 2 = 15 and R\ = R 2 = 30, reveals that translocation is 
faster when a\ > a 2 than in the reverse case when the cylinder radii are interchanged so 
that a 2 > a\. This observation may be explained by the fact that chain migration into the 
outer cylinder from the donor sphere involves a large entropic penalty, which is reduced with 
increase in the outer cylinder radius. 

The durations of the three stages of translocation are illustrated in Fig. [7J for several 
values of a 2 with a\ = 4, and vice versa, with Mi = M 2 = 15 and R\ = R 2 = 30. Figure [7] 
reveals a decrease in (n) with increase in a\ for fixed a 2 , owing to the concomitant lowering of 
the entropic barrier to translocation into the outer cylinder from the donor sphere. Varying 
a 2 with cti held fixed has no significant effect on (ti), and the slight dependence of (ti) on a 2 
arises from the restriction that the leading chain end be located at the outer cylinder-inner 
cylinder interface within a radius of a 2 from the pore axis, at the end of the first stage of 
translocation. The second stage of translocation is hastened by an increase in a 2 for fixed 
ai, owing to the reduction in the entropic cost of translocation into the inner cylinder as a 2 
is increased. Figure [7J suggests the existence of a minimum in (r 2 ) with respect to a\ for 
fixed a 2 , which occurs when at < a 2 . Varying a\ when a\ < a 2 may produce two opposing 
effects on (r 2 ). First, for a\ < a 2 , the entropic gain on entering the wider inner cylinder 
decreases with increase in a±. On the other hand, chain segments initially contained in 
the donor sphere must migrate across the outer cylinder before they can reach the inner 
cylinder, and the corresponding entropic cost is lowered as a\ is increased. These opposing 
effects may combine to produce a minimum in (r 2 ) with respect to a% for a\ < a 2 , given 
a fixed a 2 . When a\ > a 2 , the loss of entropy on entering the narrower, inner cylinder 
grows with increase in a\, leading to an increase in (t 2 ). A slight decrease in {Y3) is seen 
with increase in a\ for fixed a 2 . This decrease may be attributed to the fact that, as a\ 



16 

increases, configurations in which a greater number of chain segments is present in the outer 
cylinder at the beginning of the third stage, rather than in the donor sphere, become more 
favorable. The former configurations have a higher free energy than the latter and, hence, 
lead to a larger free energy drop upon complete chain transport into the receptor sphere, thus 
favoring translocation. A minimum in (r 3 ) with respect to a 2 for fixed a\ when a 2 < ai is 
also observed. An increase in a 2 may influence (r 3 ) in two opposing ways. As a 2 is increased, 
the entropy gained during the third stage of translocation into the receptor sphere is lowered. 
However, simultaneously, it becomes entropically more favorable for chain segments to move 
through the inner cylinder with increase in a 2 . The combination of these two effects may 
be responsible for producing a minimum in (r 3 ) with respect to a 2 . The monotonic increase 
in (r 3 ) with respect to increasing a 2 when a 2 > a\ may be attributed to the concomitant 
lowering of the entropic gain upon chain migration into the receptor sphere from the inner 
cylinder. 

Figure [8] depicts the total translocation time as a function of M 2 for Mi = 10, and vice 
versa, with a\ = 4, a 2 = 2 and R\ = R 2 = 30. An increase in M 2 at fixed Mi is accompanied 
by an increase in (r), as expected, owing to the increased cost of chain transport through 
the inner cylinder. The apparent maximum in (r) with respect to Mi for fixed M 2 and 
for Mi < M 2 may have arisen from numerical errors in our procedure (described in the 
Appendix). However, an initial decrease in (r) with increase in Mi for fixed M 2 may be 
attributed to the hastening of the second and third stages of translocation with increase in 
Mi. As Mi is increased, an increasing number of chain segments is likely to be present in the 
outer cylinder, rather than in the donor sphere, at the end of the first stage of translocation. 
Such configurations provide a greater driving force for subsequent chain transport into the 
inner cylinder and the receptor sphere, as opposed to higher-entropy configurations with a 
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large number of chain segments in the donor sphere. This argument is supported by the 
observations (not shown here) that (r 3 ) decreases with increase in M 1 or M 2 (owing to a 
concomitant increase in the entropic gain upon translocation into the receptor sphere), while 
(r 2 ) decreases as Mi is increased and increases as M 2 is increased. It should also be noted 
that an increase in (ti) occurs as M x is increased, while M 2 has no effect on (ti). The 
increase in (ti) with respect to increasing M 1 may lead to an eventual increase in (r) as Mi 
is increased at a fixed value of M 2 . 

IV. CONCLUSIONS 

In this contribution, we investigate the translocation of a Gaussian chain through com- 
posite pore geometries in the absence of excluded volume effects and hydrodynamic interac- 
tions. The translocation process is modeled as the diffusion of the translocation coordinate 
over a free energy barrier, governed by the Fokker-Planck equation subject to the radiation 
boundary condition at interfaces between pore constituents. The radiation boundary condi- 
tion allows the leading chain end to visit the interface many times before it finally migrates 
into the portion of the pore downstream of the interface. 

We illustrate our method by investigating chain translocation from one spherical chamber 
to another through a cylindrical pore and through a two-cylinder composite pore representing 
an a-hemolyin membrane channel. Our derivation of the free energy landscape is based 
upon the enumeration of all possible chain configurations during translocation, under the 
assumption that the chain length greatly exceeds the cylinder lengths. Thus, we do not 
consider situations wherein the chain completely enters the pore and subsequently diffuses 
downstream to the entrance of the receptor chamber. The latter situation may arise if the 
pore length greatly exceeds the polymer length, or if the pore width becomes comparable to 
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the polymer coil size. In such situations, the time taken by the chain to diffuse across the 
length of the pore scales with the square of the pore length. However, we do not attempt to 
quantify such effects here. Moreover, our procedure for obtaining the probability densities of 
chain configurations within the pore disallows the chain from forming hairpin configurations 
at the interface between constituent parts of the pore, and is consequently restricted to small 
pore diameters. 

Our results reveal that the time of translocation from the donor to the receptor sphere 
connected by a cylindrical pore monotonically increases with increase in the cylinder length. 
In contrast, the earlier study of Wong and Muthukumar 22 , which restricted reverse chain 
motions at the interface between the cylinder and the receptor sphere, predicted the oc- 
currence of a minimum in the translocation time with respect to cylinder length. These 
observations suggest that reverse chain motions at the interface, which are captured by the 
radiation boundary condition, may contribute significantly to the total translocation time. 
We further establish that the translocation of a chain through an a-hemolysin channel is 
faster when the chain is introduced on the cis side of the pore, rather than the trans side, as 
may be expected based on entropic considerations. This observation is consistent with the 
experimental results of Henrickson et al.— , who found the frequency of translocation-induced 
ionic-current blockades to be higher when polynucleotide chains were introduced on the cis 
side, rather than the trans side of an a-hemolysin pore. Henrickson et al. have suggested 
that pore-polymer electrostatic interactions may explain the observed asymmetry. However, 
while pore-polymer interactions may influence the polymer concentration in the vicinity of 
the pore, attractive interactions serve to greatly slow the translocation process^ 1 ^. On the 
other hand, the lower entropic cost of entering the pore on the cis side and subsequent 
chain migration may explain the observations of Henrickson et al. Finally, the approach 
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presented in this contribution may be applied to study translocation through other asym- 
metric biological or synthetic pores, including conically-shaped pores^, as well as other pore 
geometries. 



APPENDIX 

Translocation through a two-cylinder composite proceeds through intermediate configu- 
rations wherein the chain straddles the two-cylinder interface. Accordingly, we obtain the 
probability density for a chain of N segments within a cylinder of radius a and length M, 
whose one end is tethered at a distance of c from a cylinder end surface and lies within a 
radial distance of a' from the cylinder follows: 



P Cl (a, M,N,a)= > - 2 f- ^ exp - — + — — — 

M ^L x K,k Ji (^o,fc) L V a 2 M 2 J 6 



(35) 



Similarly, the probability density for a chain within the cylinder with both ends tethered at 
distances of c\ and c 2 from the cylinder end surfaces, and lying within radial distances of a' 
and a" from the cylinder axis, is 

^c 2 (a> M,J\,a ,a ) - — > — 2~ (— IJ ^ 



x exp 



Vo, k m 2 7r 2 \ N 



(36) 



a 2 M 2 J 6 

We approximate the probability density for a chain whose one end lies within 
the outer cylinder and the other end within the inner cylinder with the product 
P Cl (a 1 ,M 1 ,N 1 ,a min ) P Cl (a 2 , M 2 , N 2 , a min ) /(vra^J, where a min = min(ai,a 2 ). Because 
Eqs. (1331) and fl3"6"|) involve averaging over all possible positions of the chain segment at 
the interface between the cylinders, our results do not reduce exactly to the results of Sec. 
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HT1 upon setting a\ = a 2 . However, we have verified that our approximations qualitatively 
reproduce the behavior of the translocation time obtained from Sec. [Til for the case a\ = a 2 . 

With the aid of the above definitions, the free energies of the chain configurations illus- 
trated in Fig. [5] are defined by the following expressions: 

F ia) (I = 0) = - In [P Sl (Ri, N) Tial] (37) 

F (b) (l <Ri) = - In [P Sl (R 1} N-l) P Cl ( ai , M x , I, ai)] (38) 
F (c) (l = R l < N) = -In [P Sl (R u N-l) P C2 (a l5 M u I, a u a min )] (39) 
F {d) (l — Ri — N) — —In [P C2 (oi, Afi, N, a x , a min )] (40) 

P (e) (m < i^, i2f<JV-l) = -ln [P Sl (Ri, N — R t — m) P C2 (oi, Mi, P z , oi, a min ) 

x P Cl (a 2 , M 2 , m, a min ) / (na 2 min )] (41) 

P (/) (m = R m ,R m <N-R l ,R l <N -1) = - In [P Sl (P 1; N - Ri - R m ) 

, Oi, Omin 

) Pc 2 ( a 2, M 2 , P m , a min , a 2 ) / (TTCtmin)] ( 42 ) 

F (g)( m = R m = N - Ri, Ri < N) = - In [P C2 (a h Mi, Pj, ai, a min ) 

x Pc< 2 (a 2 , M 2 , P m , a min , a 2 ) / (vra^ in )] (43) 

F( h ){m = N - Ri) = -In [P C2 (a x , M 1; P ; , a x , a min ) P Cl (a 2 , M 2 , JV - R h a min ) / (vra^ in )] 

(44) 

F(i)(N - Ri <m < R m ) = - In [P Cl (ai, Mi, N - m, a min ) P Cl (a 2 , M 2 , m, a min ) / (7ra^ in )] 

(45) 

F(j){m = R m < N) = - In [P Cl (a u M u N - m, a min ) P C2 (a 2 , M 2 , m, a min , a 2 ) / (7ra^ in )] 

(46) 
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F {k) (m = R m = N) = - In [P Ca (a 2 , M 2 , iV, a min , a 2 )] (47) 

F (z) (l < n < N - R t - Rm) = - In [P 5l (Pi, N -R l -R m -n) P C2 (01, Mi, P z , en, a min ) 

x P C2 (a 2 , M 2 , P m , a min , a 2 ) P Sl (R 2 , n) / (vra^ in )] (48) 

F( m )(n = N -Ri- R m ) = - In [Pc 2 (a 1; Mi, P z , ai, a min ) P Ca (a 2 , M 2 , P m , a min , a 2 ) 

x P 5l (P 2 ,n)/(7ray] (49) 

F (n) (N — Ri — R m < n < N — R m ) = — In [P Cl (a x , M 1 ,N — R m — n, a min ) 

x P C2 (a 2 , M 2 , P m , a min , a 2 ) P Sl (P 2 , n) / (vra^ in )] (50) 

P (o) (n = iV - R m ) = - In [P C2 (a 2 , M 2 , P m , a min , a 2 ) P 5l (P 2 , n)\ (51) 

F w (iV - P m < n < N) = - In [P Cl (a 2 , M 2 , N — n, a 2 ) P 5l (P 2 , n)] (52) 

F (q) (n = N) = - In [P Sl (P 2 , iV) vra^] (53) 

We identify the probabilities P(P Z ) and P(P m |P z ) employed in Eqs. (P]) -(l3"4l with the 
weights of the corresponding chain configurations, as given by the expressions: 



P Sl (Pi, JV - Ri) P C2 (ai, Mi, P z , ai, a min ) , M a < P z < iV - 1 

(54) 

P C2 (ai,M l5 iV , Oi, Omin ) , Ri = N 
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and 



P Sl (Ri, N - Ri - R m ) P C2 (ai, Mi, R h ai, a n 



xP C2 (a 2 , M 2 , R m , a min , a 2 ) / (ira 2 min ) , M x < R t < N - 1, < N - R l 



P(Rm\Ri) = { 



)Pc 2 {a 2 ,M 2 ,R 

mi O-miri) 

a 2 ) I (7ra min , . 



M 1 <R l <N-l,R m = N-R l 



P Cl (a 1 ,M h N -R 

mi <^min )P C2 (a 2 ,M 2 ,R 

mi ^miri; 

a 2 ) I (na 



2 - ) 

mm/ ' 



(55) 



M l < R t < N - 1, N - Ri < R m < N or R t = N, R m < N 



Pc 2 (a 2 ,M 2 ,R 

mi C-miii) a 2 ) , M 1 <R l <N,R m = N 



The rate constants k\ and k 2 appearing in the radiation boundary conditions in Eqs. (J2 
and (|28l) may be obtained by means of the following expressions: 

- In [P S1 (R U N- R t ) P C2 (ai, Mi, R h a u a min )] 



dF (I) 
dl 



+ In [P S1 (Ri,N — Ri + 1) P Cl (ai, M x , R t - I, at)} , Ri < N 



- In [P C2 (a 1; Mi, N, a x , a min )] + In [P S1 (i?i, 1) P Cl (a 1; M t ,N — 1, <n)] 



(56) 
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- In [P Sl (RuN-Ri- R m ) P C2 (ai, M u R h a u a min ) P C2 (a 2 , M 2 , P m , a min , a 2 ) / (vra^ in )] 
+ In [P Sl (R 1 ,N-R l -R m + 1) P C2 (ai, Mi, P, , Gl, OminJ 
x P Cl (a 2 , M 2 , P m - 1, a min ) / (vra^J] ,R m < N -Ri 

~ ln [-Pc 2 Mi, a l> a min) Pc 2 («2, M 2 , P m , a min , a 2 ) / (™^ in )] 

+ ln [P 5l (R 1 ,l)P C2 (a 1 ,M 1 ,R l ) P Cl (a 2 , M 2 , P m - 1, a min ) / (vra^J] , 



- ln [P Cl (ai, Mi, N - R m , a min ) P C2 (a 2 , M 2 , P m , a min , a 2 ) / (vra^ in )] 

+ ln [P Cl (ai, Mi, iV - P m + 1, a min ) P Cl (a 2 , M 2 , P m - 1, a min ) / (7ra^ in )] , 
N — Ri < R m < N 

- ln [P C2 (a 2 , M 2 , iV, a min , a 2 )] + In [P Cl (ai, Mi, 1, a min ) 

x P Cl (a 2 , M 2 , JV — 1, a min ) / (vra^J] , R m = N 



OF (m) 



= < 



dm 



Rm — N — Ri 



(57) 



P Sl (Pi, TV - P,) P C2 (ai, Mi, P,, ai, a min ) , R t < N 



(58) 



Pc 2 (ai, Mi, TV, ai, a min ) , R t — N 
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W m {R h R m ) ~ <^ 



P 5l (Pi, N — R t — R m ) P C2 (ai, M u R h a u a min ) 



xP C2 (a 2 , M 2 , R m , a min , a 2 ) / (vra^ in ) , 



Rm < N — R, 



P C2 (ai, Mi, P z , ai, a min ) P C2 (a 2 , M 2 , P m , a min , a 2 ) / (7ra^ in ) , 



-Rm — N — Rl 



P Cl (a 1 ,M 1 ,N - R m , a min ) P C2 («2, M 2 , P m , a min , a 2 ) / (vra^J , 



N — Ri < R m < N 



02) , R m = N 



(59) 



0Wj 
9/ 



P Sl (Pi, AT - P) Pc 2 (ai, Mi, P, ai, a min ) 

-P S1 (Pi, iV - P + 1) P Cl (ai, Mi, P - 1, ai) , R t < N 
P C2 (ai, Mi, AT, ai, a min ) - P 5l (Pi, 1) P Cl (ai, M 1 , N — 1, Ql ) , Ri = N 



(60) 
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and 



dW m 
dm 



P Sl (Pi, N - Ri- R m ) P C2 (ai, Mi, P ; , ai, a min ) P C2 (a 2 , M 2 , P m , a min , a 2 ) / (vraLJ 

-P Sl (Pi, N — Ri — R m + 1) P C2 (ai, M 1 ,R h a 1 , a min ) P Cl (a 2 , M 2 , P m - 1, a min ) 
/ (ttcj^J ,R m <N-Ri 

Pc 2 (ai, M 1; P^, ai, a min ) Pc 2 (a 2 , M 2 , P m , a min , a 2 ) / (7ra^ in ) 

)P Cl (a 2 ,M 2 ,P m -l 

) / ( 7ra min) > 

-Rm = N — Ri 

P Cl {a u M u N - R m , a min ) P C2 (a 2 , M 2 , P m , a min , a 2 ) / (7ra^ in ) 

-P Cl (ai, Mi, N - R m + 1, a min ) P Cl (a 2 , M 2 , P m - 1, a min ) / (vra^ in ) , 
N — Ri < R m < N 

Pc 2 (a 2 , M 2 , AT, a min , a 2 ) - Pci (ai, Mi, 1, a min ) P Cl (a 2 , M 2 , N — 1, a min ) / (7ra^ in ) , 



P m = iV 



(61) 
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FIG. 1: Chain configurations during translocation from the donor to the receptor sphere through 
a cylindrical pore. 
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FIG. 2: (Color online) (a) Total average translocation time, (b) average time taken for the com- 
pletion of the first stage of translocation, and (c) average time taken for the completion of the 
second stage of translocation from the donor to the receptor sphere through a cylindrical pore 
as a function of cylinder length M for several values of R2 with R\ = 30 and a = 3 for a chain 
possessing ./V = 300 segments. 
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FIG. 3: (Color online) (a) Total average translocation time, (b) average time taken for the com- 
pletion of the first stage of translocation, and (c) average time taken for the completion of the 
second stage of translocation from the donor to the receptor sphere through a cylindrical pore 
as a function of cylinder length M for several values of R2 with R\ = 30 and a = 5 for a chain 
possessing N = 300 segments. 
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FIG. 4: (Color online) (a) Total average translocation time, (b) average time taken for the com- 
pletion of the first stage of translocation, and (c) average time taken for the completion of the 
second stage of translocation from the donor to the receptor sphere through a cylindrical pore as 
a function of cylinder length M for several values of a with R\ = 30 and R2 = 60 for a chain 
possessing N = 300 segments. 
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FIG. 5: Chain configurations during translocation from the donor to the receptor sphere through 
an a-hemolysin membrane channel represented as a composite two-cylinder pore. For the purpose 
of illustration, we have selected the case a\ > a.2- 
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FIG. 6: (Color online) Average time of translocation from the donor to the receptor sphere through 
an a-hemolysin channel with (a) an outer cylinder radius of a± = 3 for several values of the inner 
cylinder radius a,2 (dashed line) and vice versa (dotted-dashed line), and (b) an outer cylinder 
radius of a\ =4 for several values of the inner cylinder radius a>2 (dashed line) and vice versa 
(dotted-dashed line), and with R\ = R2 = 30 and M\ = M2 = 15, for a chain possessing N = 100 
segments. 
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FIG. 7: (Color online) Average time required for the completion of the (a) first, (b) second, and 
(c) third stages of translocation through an o-hemolysin pore having an outer cylinder radius of 
at = 4 for several values of the inner cylinder radius a<i (dashed line) and vice versa (dotted-dashed 
line), and with R\ = R2 = 30 and M\ = M2 = 15, for a chain possessing TV" = 100 segments. 
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FIG. 8: (Color online) Average time of translocation from the donor to the receptor sphere through 
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receptor sphere through a cylindrical pore as a function of cylinder length 
M for several values of a with Ri = 30 and R 2 = 60 for a chain possessing 
iV = 300 segments 



Chain configurations during translocation from the donor to the receptor 
sphere through an a-hemolysin membrane channel represented as a composite 
two-cylinder pore. For the purpose of illustration, we have selected the case 

a\ > a 2 

(Color online) Average time of translocation from the donor to the receptor 
sphere through an a-hemolysin channel with (a) an outer cylinder radius of 
ax = 3 for several values of the inner cylinder radius a 2 (dashed line) and 
vice versa (dotted-dashed line), and (b) an outer cylinder radius of a\ = 4 
for several values of the inner cylinder radius a 2 (dashed line) and vice versa 
(dotted-dashed line), and with Ri = R 2 = 30 and Mi = M 2 = 15, for a chain 

possessing iV = 100 segments 

(Color online) Average time required for the completion of the (a) first, (b) 
second, and (c) third stages of translocation through an a-hemolysin pore 
having an outer cylinder radius of a\ — 4 for several values of the inner 
cylinder radius a 2 (dashed line) and vice versa (dotted-dashed line), and with 
Ri = R 2 = 30 and Mi = M 2 = 15, for a chain possessing N = 100 segments. 
(Color online) Average time of translocation from the donor to the receptor 
sphere through an a-hemolysin channel with an outer cylinder length of Mi = 
10 for several values of the inner cylinder length M 2 (dashed line) and vice 
versa (dotted-dashed line), and with Ri = R 2 = 30, a>i = 4 and a 2 = 2 for a 
chain possessing iV = 100 segments 



